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ABSTRACT 

Context. LS 5039 has been observed with several X-ray instruments so far showing quite steady emission in the long term and no 
signatures of accretion disk. The source also presents X-ray variability at orbital timescales in flux and photon index. The system 
harbors an O-type main sequence star with moderate mass loss. At present, the link between the X-rays and the stellar wind is unclear. 
Aims. We aim here at studying the X-ray fluxes, spectra, and absorption properties of LS 5039 at apastron and periastron passages 
during an epoch of enhanced stellar mass loss, and the long term evolution of the latter in connection with the X-ray fluxes. 
Methods. New XMM-Newton observations have been performed around periastron and apastron passages in September 2005, when 
the stellar wind activity was apparently higher. April 2005 Chandra observations on LS 5039 are also revisited. Moreover, a compi- 
lation of Ho- EW data obtained since 1992, from which the stellar mass loss evolution can be approximately inferred, is carried out. 
Results. XMM-Newton observations show higher and harder emission around apastron than around periastron. No signatures of ther- 
mal emission or a reflection iron line indicating the presence of an accretion disk are found in the spectrum, and the hydrogen column 
density (N H ) is compatible with being the same in both observations and consistent with the interstellar value. 2005 Chandra obser- 
vations show a hard X-ray spectrum, and possibly high fluxes, although pileup effects preclude obtaining conclusive results. The Ho- 
EW shows yearly variations of a ~ 10%, and does not seem to be correlated with X-ray fluxes obtained at similar phases, unlike it 
would be expected in the wind accretion scenario. 

Conclusions. 2005 XMM-Newton and Chandra observations are consistent with 2003 RXTEfPCA results, namely moderate flux and 
spectral variability at different orbital phases. The constancy of the Njj seems to imply that either the X-ray emitter is located at 
^ 10 12 cm from the compact object, or the real density in the system is 3 to 27 times smaller than the one predicted by a spherical 
symmetric wind model. We suggest that the multiwavelength non-thermal emission of LS 5039 is related to the observed extended 
radio jets and unlikely produced inside the binary system. 

Key words. X-rays: binaries - stars: individual: LS 5039 - Radiation mechanisms: non-thermal 
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* The Ha results presented here are based on observations made 
with ESO Telescopes at the La Silla or Paranal Observatories un- 
der program IDs 67.D-0229(A), 69.D-0628(A) and 075.D-0591(A); 
the Observatoire de Haute-Provence; the Observatorio do Pico dos 
Dias / LNA, Brazil; the G.D. Cassini telescope operated at the Loiano 
Observatory by the Osservatorio Astronomico di Bologna; and the 
Nordic Optical Telescope, operated on the island of La Palma jointly 
by Denmark, Finland, Iceland, Norway, and Sweden, in the Spanish 
Observatorio del Roque de los Muchachos of the Instituto de Astrofisica 
de Canarias (those data were taken using ALFOSC, which is owned 
by the Instituto de Astrofisica de Andalucfa (IAA) and operated at 
the Nordic Optical Telescope under agreement between IAA and the 
NBIfAFG of the Astronomical Observatory of Copenhagen). 



1. Introduction 

LS 5039 (V479 Set) is an X-ray binar y (Mo tch et al. [1997b 
located at 2.5+0.5 kpc (Casares et al. 120051 C05 hereafter). 
The orbital ephemeris of the source have been obtained by 
C05, being: T =HJD 2451943.09 + 0.1, periastron passage 
-phase 0.0-; P orb = 3.90603 ± 0.00017 days, orbital period; 
e = 0.35 + 0.04, eccentricity; and i ~ 15°-60°, inclination^ 
LS 5039 presents — 1—100 milliarcseconds- radio jets (Paredes 
et al. 120001 120021) . shows X-ray variability on timescales sim- 
ilar to the orbital one (Bosch-Ramon et al. 120051 BR05 here- 
after), and has been detected at very high-energy gamma-rays 
(Aharonian et al. 120051) . which virtually confirms its association 
with a CGROfEGKET (EGRET from now on) source (Paredes 
et al. 120001 ). Interestingly, the TeV emission varies with the or- 



1 We note that all through this work, the errors in phases -not shown- 
range between 0.03-0.04 due to the uncertainties in To and P or t,. 
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bital period (Aharonian et al. 2006). Recently, LS 5039 has been 
also detected by INTEGRAL (Goldoni et al. |2DD6] l up to 100 keV 
(De Rosa et al. 2006 ). The nature of the compact object is still 
uncertain. C05 suggest that it may be a black hole, although there 
is an on-going debate on this issue, and some authors have pro- 
posed that LS 5039 is in fact a young non-accreting pulsar (see, 
e.g., Martocchia et al. |2005l -M05 hereafter-; Dubus 12006] ) . 

All X-ray observations performed so far by imaging instru- 
ments failed to show strong flux and photon index variations. 
RXTE/PCA (RXTE from now on) observed also the source in 
1998 (Ribo et al. [19991 and 2003 (BR05), although the latter au- 
thors found that RXTE data were likely contaminated by diffuse 
X-rays from the Galactic Ridge. Still, 2003 RXTE data showed 
flux variations by a factor of ~ 2 unrelated with diffuse back- 
ground emission, peaking smoothly around phase 0.8 and more 
sharply at other phases. These peaks were apparently accompa- 
nied by spectral hardening. At phases ~ 0.8, higher and harder 
emission has been also observed at TeV energies (Aharonian 
et al. 120061 ). and simultaneous Chandra observations have ap- 
parently shown a similar behavior in X-rays (Horns et al. 2006). 
A possible relationship between the mass-loss rate of the com- 
panion star and the X-ray emission in LS 5039 was discussed 
in previous works and a correlation suggested in the context of 
wind accretion powering the X-rays (e.g., Reig et al. 2003 -R03 
hereafter-; McSwain et al. l2004l -MS04 hereafter-). 

In this work, we present recent XMM-Newton observations 
of LS 5039 carried out in September 2005 during periastron and 
apastron passages, in an epoch of inferred high stellar mass-loss 
rate. Chandra data taken in the same year is reanalyzed. In ad- 
dition, we investigate the long-term Ha EW evolution, linked to 
the evolution of the stellar mass-loss rate, and compare it with 
the X-ray fluxes looking for possible correlations. These results 
are put in context with information at other wavelengths and a 
physical scenario is suggested. 

2. X-ray observations 

2.1. XMM-Newton observations 

The XMM-Newton observations presented here were performed 
in two runs, on the 22th and 24th of September 2005, with EPIC 
pn on times of 15.8 ks and 10.4 ks (very similar to those of 
MOS1 and MOS2), respectively (see Table [TJ. The correspond- 
ing observed phases are 0.49-0.53, around apastron, and 0.02- 
0.05, right after periastron. These observations were triggered in 
the context of a target of opportunity (ToO) program linked to 
the strength of the Ha line, which gives information on the state 
of the stellar wind. The Ha EW was found in late August 2005 
to be above of -2.4A, defined as the threshold set to trigger the 
XMM-Newton observations and significantly larger than the val- 
ues observed during the previous years. XMM-Newton observa- 
tions were carried out in September due to visibility constraints. 
The physical motivation of the ToO was that, if direct accretion 
of the stellar wind is taking place in the system, the detection 
of higher X-ray fluxes, and/or accretion spectral features like a 
multi-color black-body and/or a reflection iron line, seems more 
likely during epochs of larger stellar mass loss, i.e. higher accre- 
tion rate. The convention for the Ha EW in this work is to take 
this value as negative when in absorption, i.e. it increases in the 
same sense as the stellar mass-loss rate. 

We focus here on the data obtained with the EPIC instru- 
ment, since the X-ray spectrometer (RGS) data did not have 
enough statistics due to the faintness of the source. The MOS 
and pn cameras were in small and full window respectively, with 



Table 1. X-ray observations and spectral results. 





Date 


2005-04-13 


2005-09-22 


2005-09-24 


Mission 


Chandra 


XMM-Newton 


XMM-Newton 


OBSID 


6259 


0202950201 


0202950301 


Orbital phase 


0.86-0.87 


0.49-0.53 


0.02-0.05 


MJD start 


53473.116 


53635.714 


53637.779 


MJD stop 


53473.174 


53635.897 


53637.899 


/, {L x ) 


2±1 (1.5) 


1.18!°$ (0.9) 


0.74±0.07 (0.6) 


Photon index 


1.1±0.2 


1.51 ±0.02 


1.59±0.03 


Nh (phabs) 


6±1 


6.3±0.1 


6.2±0.2 



Notes: fluxes and luminosities are given unabsorbed in the 0.3-10 
keV energy band in units of 10~" erg cirr 2 s -1 and 10 34 erg s -1 , 
respectively. The Nh units are 10 21 cirr 2 . Quoted errors are at the 
1 cr confidence level. Ephemeris taken from C05. 



medium filters. The reduction and analysis of the data have been 
performed following the standard procedure using the SAS 7.0.0 
software package to filter data and create spectra and lightcurves, 
and Xspec 12.3.0 for the spectral analysis. No large background 
flare was present in any of the two data sets. The 0.3-12 keV 
count rates per camera around phases 0.5 and 0.0 were 0.46 and 
0.29 cts s 1 for MOS1, 0.47 and 0.30 cts s 1 for MOS2, and 1.09 
and 0.71 cts s _1 forpn, respectively. 

The lightcurves for the two runs, adding the data from pn, 
MOS1 and MOS2 cameras and taking time bins of 500 s, are 
presented in Fig. Q] for periastron (bottom) and apastron (top) 
passages, respectively. Around apastron, the emission is variable 
on timescales of a few hours. A constant fit to the data gives a 
reduced^ squared Cfcd d ; 27 degrees of freedom -d.o.f.-) of 6.8, 
and the maximum count rate variation is a 25% during this run, 
with deviations of « ±5<x from the mean value. Around peri- 
astron, the lightcurve is still somewhat inconsistent with being 
constant (xL d = 2.4; 18 d.o.f.), with a maximum count rate vari- 
ation of 25% and deviations of « ±2<x from the mean value. 
We studied also the hardness ratio evolution for both runs, which 
is compatible with being constant. Searching for pulsations, no 
significant period peaks are found in the frequency range of 0.01 
to 83 Hz and the pulsed fraction has been constrained to be less 
than =s 10%, similar to the results reported by M05 for 2003 
XMM-Newton data. 

We show in Fig. [2] the spectrum of the source around phase 
0.5, altogether with the deviations in cr of the best fit model 
from the data. We have added the data obtained from the pn, 
MOS1 and MOS2 cameras. The spectrum around phase 0.0 is 
not shown though it looks very similar to the one around phase 
0.5. The spectrum of the emission is well represented by an ab- 
sorbed power-law Cf r 2 ed = 1.08 -637 d.o.f.- and 0.93 -312 d.o.f. - 
for phases 0.5 and 0.0, respectively). The unabsorbed fluxes, 
photon indices (0.3-10 keV) and A^h values for the two runs 
are given in Table Q] The flux between periastron and apastron 
increases by a 60% at a level of 4.4cr, and the photon index 
rises up marginally by 0.08 (2.2<x). This behavior, despite of be- 
ing marginal, give support to the RXTE finding that the higher 
the flux, the harder the spectrum. For the absorption model, we 
used phabs, which assumes updated photoelectric cross-sections 
from Balucinska-Church & Mc Cammon ( 1992 ). The A^h values 
are compatible within errors between both runs. 

A multi-color black-body (bbody) and a Gaussian line com- 
ponent have been added, separately, to the absorbed power-law 
model for the two runs. We obtain upper-limits for the multi- 
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Fig. 1. Lightcurves at periastron (top) and apastron (bottom). 
The average value is given as a long-dashed line. We note that 
the time axis range differently in each lightcurve. 




Fig. 2. Spectral data from pn (upper spectrum), MOS1 and 
MOS2 (lower spectra) cameras, and the best fit to the three of 
them: an absorbed power-law, around phase 0.5 (those around 
phase 0.0 look very similar and are not shown here). The devi- 
ations in <x of the model from the data are shown in the lower 
panel. 



color black-body flux in the two runs similar to those found by 
M05. At the same phases, we look for flux upper-limits of a pos- 
sible iron line and obtain an upper-limit for the flux of about a 
2% the total one. The lack of a multi-color black-body and a line 
component in the data is consistent with previous imaging in- 
strument observations and confirms that the iron line reported in 
Ribo et al. ( 1999) was due to background contamination. 

2.2. Chandra data 

In April 2005, Chandra observed LS 5039 in the ACIS-S mode 
(time frame: 1.6 s) during 5 ks (see Table [TJ at a time corre- 
sponding to the orbital phase 0.86-0.87, at which RXTE ob- 
served a smooth maximum (BR05). The Chandra Interactive 
Analysis of Observations software package (CIAO 3.3.0.1) has 
been used to extract the spectrum and the lightcurve, following 
the standard procedure given in the Chandra analysis threads. 



Sherpa (CIAO 3.3) has been used to analyze the spectrum. The 
ACIS count rate was relatively high, 0.3 cts s (0.3-10 keV). 
Chandra data is likely moderately affected by pileup (see The 
Chandra Proposers' Observatory Guide, Fig. 6.19). Therefore, 
we have used the pileup model available in Sherpa to obtain our 
fits (Davis[200l). 

The spectrum can be well fitted by a power-law model with 
absorption {xsphabs, the sherpa version of phabs) and pileup. 
We perform the spectral fit fixing the event pileup fraction pa- 
rameter to 0.95 to avoid unreasonably high values for the flux 
(i.e. ~ 10~ 10 erg cirT 2 s _1 ). We obtain a relatively low grade mi- 
gration parameter of 0.2 + 0.1. The pileup fraction is « 11%. 
Therefore, the effect of pileup on these observations seems to be 
only moderate. The results of our fit are given in Table Q] The 
photon index value, harder than that found in almost all previous 
observations, has been also calculated without a pileup model 
by removing the central pixels, i.e. the most affected by pileup 
(pileup effects can indeed produce an artificial hardening of the 
spectrum), obtaining compatible spectral results. No disk or line 
components are required to fit the data. The lightcurve is consis- 
tent with being constant. 



3. Ha line observations 

We present in Table 2 the list of the journal of optical observa- 
tions and the corresponding Ha EW. Most of these optical spec- 
troscopic observations of LS 5039 were carried out by us at reg- 
ular intervals of time between September 2000 and July 2002, 
with the aim to construct an orbital solution and to monitor the 
variations of the EW of the Ha line, the most convenient line 
to study the wind of LS 5039. A new observing campaign was 
initiated in 2004 in order to determine the time at which the Ha 
EW would go above -2.4 A, the condition required to trigger the 
2005 XMM-Newton observations. 

Various instruments were used in the course of this extended 
campaign. At the Observatoire de Haute-Provence (OHP), 
we observed most of the time with the ELODIE Echelle 
Spectrometer. Some additional spectra were obtained with the 
CARELEC spectrograph (Lemaitre et al. 1990). At ESO, the 
majority of the data were acquired with the FEROS0 Echelle 
Spectrometer (Kaufer et al. 1999). FORS1 observed the source 
on two occasions in 2001. LS 5039 was also monitored over 
an orbital cycle with the ESO-NTT and EMMI during a dedi- 
cated program in July 2002. Finally, the spectrum obtained in 
1992 is described in Motch et al. (1997). The spectra obtained 
on 1999-08-18 and 2004-07-19 were acquired at Osservatorio 
Astronomico di Bologna with the 1.52 m Loiano telescope 
and the BFOSC instrument and at the NOT we used the 
ALFOSC spectrograph. At the Brazilian National Observatory 
of the Observatorio do Pico dos Dias, we used the ESPCASS 
Cassegrain grating spectrometer mounted on either the 0.6 m or 
the 1.6 m telescope. Instrument settings are summarized in the 
Table 2. 

The Ha line has an intrinsic FWHM of ~8 A and is therefore 
relatively well resolved in most of our observations. Due to the 
southern location of the object, observations from the OHP and 
Canary Islands were always acquired at a high air mass, typi- 
cally 1.9 at OHP. At these high air masses, several telluric ab- 
sorption lines start to be noticeable and can substantially change 
the value of the Ha EW, either by directly adding to the absorp- 
tion profile or by changing the value of the reference continuum. 



2 All FEROS observations were obtained during the Brazilian time 
allocation. 
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Table 2. Ha EW measurements obtained from 1992 to 2006*. 



Date 

(yyyy-MM-dd) 


MJD 


Orb. phase 


Observatory 


Instrument 


Resolution 


Setting** 


Ho- EW 
(A) 


Error 
(A) 


1992-04-19 


48731.250 


0.85 


ESO-2.20 


EFOSC2 


3,000 


- 


-2.80 


0.07 


1993-06-19 


49257.010 


0.45 


OHP-193 


CARELEC 


450* 




-2.50 


0.10 


1999-08-18 


51408.842 


0.35 


OAB 


BFOSC 


1,600 


Grism#7 


-2.60 


0.10 


2000-09-08 


51795.813 


0.42 


OHP-193 


ELODIE 


45,000 




-3.24 


0.10 


2000-09-10 


51797.768 


0.92 


OHP-193 


ELODIE 


45,000 


- 


-2.91 


0.15 


2000-09-12 


51799.768 


0.43 


OHP-193 


ELODIE 


45,000 




-3.23 


0.15 


2000-10-18 


51835.984 


0.70 


ESO-1.52 


FEROS 


48,000 


- 


-2.64 


0.05 


2000-10-19 


51836.983 


0.96 


ESO-1.52 


FEROS 


48,000 


- 


-2.61 


0.05 


2000-10-20 


51838.000 


0.22 


ESO-1.52 


FEROS 


48,000 


- 


-2.71 


0.05 


2001-05-24 


52053.424 


0.37 


ESO-VLT-U1 


FORS1 


990 


GRIS.600V 


-2.69 


0.07 


2001-05-26 


52055.417 


0.88 


ESO-VLT-U1 


FORS1 


990 


GRIS.600V 


-2.69 


0.07 


2001-08-25 


52146.858 


0.29 


OHP-193 


ELODIE 


45,000 




-2.95 


0.10 


2001-08-26 


52147.857 


0.55 


OHP-193 


ELODIE 


45,000 


- 


-2.98 


0.10 


2001-08-28 


52149.848 


0.06 


OHP-193 


ELODIE 


45,000 


- 


-2.91 


0.20 


2001-08-30 


52151.822 


0.56 


OHP-193 


ELODIE 


45,000 


- 


-2.92 


0.20 


2001-10-06 


52188.978 


0.08 


ESO-1.52 


FEROS 


48,000 


- 


-2.64 


0.10 


2001-10-07 


52189.992 


0.33 


ESO-1.52 


FEROS 


48,000 




-2.79 


0.10 


2001-10-08 


52191.050 


0.60 


ESO-1.52 


FEROS 


48,000 


- 


-2.50 


0.15 


2001-10-09 


52192.000 


0.85 


ESO-1.52 


FEROS 


48,000 




-2.53 


0.10 


2002-06-14 


52439.340 


0.17 


ESO-NTT 


EMMI 


5,000 


- 


-2.81 


0.05 


2002-07-08 


52463.335 


0.31 


ESO-NTT 


EMMI 


5,000 


- 


-2.78 


0.05 


2002-07-09 


52464.141 


0.52 


ESO-NTT 


EMMI 


5,000 


- 


-2.80 


0.05 


2002-07-10 


52465.052 


0.75 


ESO-NTT 


EMMI 


5,000 




-2.72 


0.10 


2002-07-11 


52466.214 


0.06 


ESO-NTT 


EMMI 


5,000 


- 


-2.70 


0.05 


2004-07-13 


53199.923 


0.89 


OHP-193 


CARELEC 


900* 


Grating 133/375 


-2.60 


0.10 


2004-07-14 


53200.937 


0.15 


OHP-193 


CARELEC 


900* 


Grating 133/375 


-2.54 


0.10 


2004-07-16 


53202.936 


0.66 


OHP-193 


CARELEC 


900* 


Grating 133/375 


-2.44 


0.10 


2004-07-18 


53204.906 


0.17 


OAB 


BFOSC 


2,000 


Grism#8 


-2.47 


0.10 


2004-09-02 


53250.935 


0.95 


LNA-0.60 


ESPCASS 


1,460 


- 


-2.60 


0.10 


2005-04-07 


53467.299 


0.34 


LNA-1.60 


ESPCASS 


3,600 


Grating 1200/750 


-2.83 


0.10 


2005-04-07 


53467.299 


0.34 


LNA-1.60 


ESPCASS 


3,600 


Grating 1200/750 


-2.77 


0.10 


2005-07-31 


53583.748 


0.16 


ESO-NTT 


EMMI 


1,100 


Grism#5 


-2.55 


0.05 


2005-08-20 


53602.935 


0.07 


LNA-1.60 


ESPCASS 


2,500 


Grating 900/550 


-2.39 


0.05 


2005-08-20 


53602.963 


0.08 


LNA-1.60 


ESPCASS 


2,500 


Grating 900/550 


-2.34 


0.05 


2005-08-21 


53603.078 


0.10 


LNA-1.60 


ESPCASS 


2,500 


Grating 900/550 


-2.43 


0.05 


2005-09-30 


53643.005 


0.33 


LNA-1.60 


ESPCASS 


730 


Grating 300/640 


-2.51 


0.10 


2006-05-07 


53862.181 


0.44 


NOT 


ALFOSC 


3,000 


Grism#16 


-2.40 


0.10 



* Spectra with EW values corrected for the presence of blended telluric lines. 
** The setting of the instrument is shown when relevant. 



These lines are clearly seen in the high resolution ELODIE in- 
dividual spectra. In the AA range 6540-6580 A, the most impor- 
tant telluric lines are located at AA 6543.89, 6546.61, 6552.62, 
6557.16, 6564.19, 6572.06 and 6574.83 A. Guided by the atmo- 
spheric transmission lines atlas of Hinkle et al. (2003 ), we could 
easily remove the sharp lines due to telluric absorption from the 
high resolution ELODIE spectra and thus fully correct for at- 
mospheric effects. Some of the FEROS spectra also exhibited 
evidences of well resolved atmospheric lines at a much lower 
level than at OHP and could be removed as well. 

At a resolution of ~ 1000 or less, however, atmospheric lines 
are heavily smeared and cannot be simply removed. They appear 
as shoulders in the blue and red wings of the broad Ha absorp- 
tion line of LS 5039 or as slight distortion in the central line 
profile. We estimated the effect of these lines on the measured 
equivalent widths by scaling the absorption spectrum of Hinkle 
et al. ( 120031 ) so that they could match the equivalent widths of 
the narrow lines observed in the mean OHP ELODIE spectrum 
and applying it to a Gaussian Ha line profile fitting that of LS 
5039. Because of the low altitude and relatively high humidity 



level, the atmospheric effects at OHP are noticeably large. We 
find that telluric absorptions can easily add 0.4 A to the equiva- 
lent width of a low resolution spectrum at airmass 1.9. Low res- 
olution spectra obtained at OHP were thus rectified using con- 
tinuum regions clear of known atmospheric lines and their EW 
were then corrected by 0.4 A. In the absence of clear measure- 
ments or estimates of water vapor column we did not try to cor- 
rect the other measurements obtained at lower air masses and 
under much dryer conditions. 

4. Discussion 

4.1. Evolution of the Ha EW and the X-rays 

We have studied the available Ha data (this work; C05; MS04; 
R03) and found, as already suggested in previous works (e.g. 
McSwain & Gies |2002l R03), that the Ha EW is not compatible 
with being constant (xLa - 3-7; 96 d.o.f.), with a weighted mean 
value of -2.66 + 0.01 A, where the error has been computed us- 
ing the error propagation method. All errors quoted here are 1 cr. 
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The weighted standard deviation of the sample is 0.21 A, giving 
an idea of the dispersion degree of the individual points. The 
yearly averages with their corresponding standard deviations, 
altogether with the individual points, are shown in Fig. [3] For 
those years when only one measurement is available, we show 
just the individual point with its error. The yearly mean errors are 
typically 0.02 -0.1 (not shown in the plot), being significantly 
smaller than the yearly standard deviations (shown in the plot), 
0.04 - 0.29 A. This relatively large dispersion seems to indicate 
that the Ha EW also varies at scales shorter than one year. Due to 
the heterogenous origin of the data, we have done a further test 
to check whether the Ho- EW is indeed variable using only data 
given in the Table 2, treated homogeneously, excluding those 
points that may be affected by atmospheric lines, concretely 
those obtained with CARELEC at OHP, and the two much ear- 
lier 1992 and 1993 points. The resulting mean is -2.62 + 0.01 A 
with a standard deviation of 0.26, being also this data subset in- 
compatible with a constant (xL d = 5.6, 24 d.o.f.). Although the 
strong dispersion of the data prevents to make strong claims, it is 
worthy noting that the Ha EW seems to increase on average dur- 
ing the period 2003-2006, which would be evidence of a higher 
stellar mass-loss rate. 

As noted in Sect. Q] a correlation between the Ha EW and 
the X-ray flux has been proposed in the past (e.g. R03; MS04) 
suggesting that larger stellar mass loss rates are correlated with 
episodes of higher X-ray luminosity. These authors included in 
their comparison non imaging 1998 RXTE data, and the used 
X-ray observations were taken at different orbital phases, being 
likely affected by additional orbital variability. Trying to avoid 
orbital variation effects and excluding RXTE data, we use XMM- 
Newton data taken in 2003 and 2005 around apastron and find 
that the flux in 2005 is only a « 10% larger (2cr level) than that 
observed in 2003. Otherwise, the Ho- EW appears larger by a 
a 14% (6tr level; between -2.76 + 0.05 A to -2.43 + 0.02 Afl 
which could indicate variations of the stellar mass-loss rate by 
a factor of ~ 2 (see, e.g., Puis et al. 119961 Table 7). Therefore, 
similar phase/epoch X-ray/Ha EW data point to a quite constant 
X-ray flux versus a significantly changing wind state. This ab- 
sence of clear X-ray flux response to a significant change in the 
stellar mass loss rate is at variance with the relation suggested by 
R03 and MS04, although it is presently just a two-points com- 
parison. 

Finally, we show in Table [3] the 0.3-10 keV unabsorbed 
fluxes, photon indices and A^h (using the wabs model) values ob- 
tained by all instruments with imaging capabilities and working 
in the mentioned energy range. Using the values provided in that 
table, we obtain the following weighted mean values for the 0.3- 
10 keV flux, (1.05 + 0.02) x 10" 11 erg cnr 2 s" 1 (fit-to-constant 
Xred = 8 - 4 ' 7 dof X the photon index, 1.54 + 0.01 Cr r 2 ed = 2.7; 
7 d.o.f.), and N H , (6.8 + 0.1) x 10 21 cm" 2 (x 2 Kd = 0.7; 7 d.o.f.). 
Therefore, the flux found by all the instruments with imaging ca- 
pabilities is clearly variable (dominated by orbital variability?) 
even when excluding the 2005 Chandra data point, the photon 
index seem to change as well with time, and the A^h is clearly 
compatible with being constant. 



4.2. The N H value and the location of the X-ray emitter 

Photo-electric absorption in the dense wind of the primary in 
LS 5039 should increase the observed Nh to values above that 
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Fig. 3. The year averages (black circles) and the individuals 
points (white squares) for the Ha EW. For those years with sev- 
eral data points, the year average points are shown and their error 
bars represent the dispersion of values. The shown data come 
from the Table 2 and previous works (C05: year 2003; MS04: 
years 1998, 1999, 2000, 2003; R03: years 1999, 2002). 



due to the interstellar medium (ISM) alone. This effect could be 
detectable when comparing the values of the A^h inferred from 
X-ray and optical observations, respectively. However, the Nh 
values inferred from XMM-Newton data (using either wabs or 
phabs; see Sect. |2]l are similar to, and somewhat smaller than, 
the value of the ISMA^h alone, which is (7.3 + l.l)x 10 21 crrT 2 , 
derived using the £ B -v value from MS04 and the ISM A^h- 
£b-v relationship given by Predehl & Schmitt ( [T995T T. Taking 
into account the optical and XMM-Newton Nh values and er- 
ror bars, we derive lower- and upper-limits of » 4.1 x 10 21 and 
6.7xl0 21 cm~ 2 for the interstellar and total A^h, respectively. This 
yields an upper-limit for the intrinsic A^h of « 2.6 x 10 21 cm -2 . 
In addition, the XMM-Newton Nh values around phases 0.5 and 
0.0 are compatible with being constant, unlike what may be ex- 
pected given the orbital eccentricity and inclination of the system 
(see Sect.Q}, which should introduce orbital variations in the A^h 
intrinsic to the wind. The 3cr-error of the difference of A^h be- 
tween apastron and periastron is =s 10 21 cm 2 , which we adopt 
as an upper-limit to the intrinsic A^h variation between these two 
orbital phases. 

To quantify whether the contribution of the stellar wind to 
the total A^h depending on the X-ray source location, we have 
modeled the wind as spherically symmetric. We have neglected 
X-ray irradiation effects since the modest X-ray luminosity of 
the source (Z> 10 34 erg s _1 ) is unlikely to ionize a large fraction 
of the wind, which could reduce the intrinsic A^h value. We es- 
timate the distance at which X-rays would strongly ionize the 
wind to be ^ 1/10 the orbital one (see, e.g., Blondin 1994). 
MS04 obtained similar constraints on the ionized wind region 
size. The radial wind velocity at a distance r from the star has 
been assumed to follow the law: V w = Voo(l - R*/rf, where 
Vca - 2440 km s is the radial velocity of the wind in the in- 
finity, /3 = 0.8 is the wind profile exponent, R, = 9.3 R Q is 
the stellar radius and r is the distance to the star (MS04; C05). 
The stellar mass-loss rate is taken to be in the range M w ~ 3.7— 
7.5 x 10~ 7 M yr~' (C05), which allows us to compute the wind 
density: p w = M w /4^r 2 V w . 

Integrating the wind density from the compact object loca- 
tion along the observer line of sight up to distances where the 
wind influence is negligible, we obtain an intrinsic A^h in the 



3 These has been obtained averaging few months data around the 
XMM-Newton observation dates, which implies to assume a steady Ho- 
behavior at these timescales. 



4 To obtain a conservative error for the Nh - Ay relationship given in 
Predehl & Schmitt ( 19951, we use here the dispersion and not the error 
of the mean. 
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Table 3. Summary of the results for all the X-ray observations on LS 5039 performed by instruments with imaging capabilities. 
The unabsorbed fluxes in the energy range 0.3-10 keV, photon indices, and A^h (with wabs model) values are given. BeppoSAX flux 
errors are the flux dispersion accounting for the count rate changes by a factor of two during the observation and neglecting photon 
index variations. ASCA flux errors are assumed to be a 10% of the flux. The given 2005 XMM-Newton Nn values are computed with 
the wabs model for comparison with 2003 values. For more details, see M05 (2003 XMM-Newton data), R03 {BeppoSAX data), and 
Yamaoka, K., private communication (ASCA data). 



Date 

yyyy-MM-dd 


MJD 


Mission 


Phase range 


Flux (0.3-10 keV) 
xl0~" erg crrr 2 s -1 


Photon index 


Nn (wabs) 
xlO 21 cm- 2 


1999-10-04 


51455.9 


ASCA 


0.38-0.47 


0.80±0.08 


1.6+0.1 


7±1 


2000-10-01 


51818.5 


BeppoSAX 


0.97-0.21 


0.5±0.2 


1.8±0.2 


io + J 4 


2002-09-10 


52527.3 


Chandra 


0.71-0.73 


0.8±0.3 


1.1+0.2 


6±2 


2003-03-08 


52706.3 


XMM-Newton 


0.54-0.57 


1.03±0.07 


1 56 +0 - 02 


7 9+0.3 
' • Z, -Q,5 


2003-03-27 


52725.8 


XMM-Newton 


0.55-0.58 


0.97±0.07 


i 4Q+0.05 


6 9 +0 - 5 

-0.3 


2005-04-13 


53473.1 


Chandra 


0.86-0.87 


2+1 


1.1+0.2 


6±1 


2005-09-22 


53635.8 


XMM-Newton 


0.49-0.53 


i i O+0.03 


1.51+0.02 


6.7±0.2 


2005-09-24 


53637.8 


XMM-Newton 


0.02-0.05 


0.74±0.07 


1.59±0.03 


6.6±0.2 



ranges ~ 8-29 x 10 21 and ~ 2-5 x 10 21 cm -2 around phases 0.0 
an 0.5, respectively. To look for these minimum and maximum 
values of the wind A^h, we have adopted the most extreme val- 
ues for the inclination and M w . The wide range of possible val- 
ues comes from the factor of 2 difference in the mass-loss rate 
and from the two extreme inclination angles. The A^h difference 
between the two orbital locations is expected to be in the range 
~ 3-27 x 10 21 cirT 2 , taking for this comparison the two extreme 
M„ values, due to possible orbital variability, but the same i (the 
one rendering either the smallest or the largest A^h difference). 
Now, the wide range of Nh values is mainly explained by the 
fact that the line of sight can cross very different density regions 
at apastron and periastron passage depending on i. 

The observed A^h seems fairly low taking into account the 
absorption from the ISM plus the intrinsic component computed 
by us. At periastron passage, the intrinsic expected A^h value is 
3 to 12 times larger than the stated upper-limit for this com- 
ponent. At apastron, the expected intrinsic A^h value is similar 
to its upper-limit. Stronger disparities are found when compar- 
ing the calculated intrinsic A^h variation between phases 0.5 and 
0.0: a factor between 3 to 27 larger than the upper-limit. It is 
worth noting that a similar wind model seems to work when ap- 
plied, for instance, to the orbital A^h variation in 4U 1538-52 
(Mukherjee et al. [2007), an X-ray binary system that presents 
X-ray absorption conditions along the orbit similar to those of 
LS 5039. Nevertheless, there are some cases in the literature 
about winds in massive isolated stars that appear to be less dense 
than expected (e.g. Kramer et al. 2003, Cohen et al. 2006]). In 
these cases, wind clumping has been invoked to explain the low 
densities of the hydrogen column density inferred from X-ray 
observations, although a very special wind structure would be 
required to prevent the observed constancy of the A^h for such 
different orbital phases in LS 5039 if the emitter is well inside 
the system. Such special wind structure may be produced by the 
slow compact object motion through the wind, although it can 
hardly have strong impact since the supersonic wind will effi- 
ciently convect outwards any possible perturbation. Another rea- 
son for this low intrinsic and/or constant A^h could be that X-rays 
are emitted where the stellar wind is already diluted. In particu- 
lar, an emitter located along the line of sight should be beyond 
one orbital radius, »2x 10 12 cm, away from the compact ob- 
ject to explain the constancy of the A^h or its low intrinsic value. 
Different orientations of the system/emitter/observer geometry 
give similar constraints, thus we adopt ~ 10 12 cm as a lower- 



limit for the emitter-compact object distance. Finally, to reduce 
the inferred A^h, the X-ray emitter may be as large as the whole 
system, although in such a case the stellar wind would probably 
be affected at much larger extent than it seems to be (e.g. via 
ionization). 



4.3. The non-thermal emitter 

The unabsorbed X-ray spectrum of LS 5039 is well explained 
by a pure power-law pointing to a non-thermal origin, and the 
data at higher energies obtained by INTEGRAL (De Rosa et al. 
2006) show no hints of cutoff up to 100 keV, unlike it would 
be expected in a thermal comptonization scenario. A location of 
the X-ray source at 10 12 cm from the compact object would 
also support a non thermal nature for the high energy emission. 
Plausible mechanisms that could produce the keV-MeV spec- 
trum could be synchrotron or inverse Compton processes, al- 
though solving this question is out of the scope of this work. 
We note that the lack of thermal features in the X-ray spectrum 
has been the main argument against accretion as the powering 
source in LS 5039, and a young non-accreting pulsar, producing 
the high energy emission via stellar/pulsar wind collision, has 
been proposed as the non thermal emitter (see, e.g., M05; Dubus 
120061 

The emission at TeV energies in LS 5039 is unlikely pro- 
duced well inside the system since photon-photon absorption 
does not seem to be the only effect modulating the TeV emis- 
sion (Aharonian et al. 12006) . In addition, the optically thin na- 
ture of the radio emission at least down to ~ 1 GHz (Marti 
et al. 1 1 9981 ). and the relatively small radio variability of the core 
(^ 50%), preclude a too compact radio emitter. We note that all 
this fits well in a scenario in which the X-ray emitter is at some 
distance from the compact object. Extended radio emission has 
been imaged in the range ~ 1-100 mas, showing an approxi- 
mately steady orientation of + 15° at different epochs and spatial 
scales as well as strong collimation (Paredes et al. 2000, 2002). 
It seems reasonable to associate such extended radio structure 
with the multi wavelength emitter. Whether this jet-like structure 
is powered by a different accretion mechanism, whatever its ra- 
diative properties, or it can be explained in the context of a more 
general young non-accreting pulsar scenario, remains an open 
question. In any case, understanding the complex radio, X-ray 
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and TeV spectra and lightcurves, and the unclear X-rays-stellar 
wind link, requires further and accurate data plus modeling. 

5. Summary 

The 2005 XMM-Newton observations of LS 5039 presented here 
were carried out around apastron and periastron passages dur- 
ing an epoch of enhanced stellar mass loss rate. The fluxes, 
photon indices, and N H are (1.18+°°^) and (0.74 ± 0.07) x 
10- 11 erg cm" 2 s -1 , 1.51 ±0.02 and L59±0.03, and (6.3 ± 0.1) 
and (6.2 ± 0.2) x 10 21 crrT 2 , at orbital phases a 0.5 and 0.0, re- 
spectively. The reanalysis of 2005 Chandra observations yields 
a flux of (2 ± 1) x 10~ n erg crrT 2 s , a photon index of 1.1+0.2, 
and N H = (6.2 ± 0.2) x 10 21 cm" 2 . The XMM-Newton fluxes 
and photon indices show variability (more marginal for the lat- 
ter), whereas the Nh is compatible with being constant. Chandra 
data point to a hard spectrum, and possible higher fluxes, al- 
though the impact of pileup adds uncertainty to these results. 
Collected Ha EW data from 1992 to nowadays show that this 
quantity varies on timescales longer, and perhaps also shorter, 
than one year. Comparison of the Ha EW at the epochs when 
2003 and 2005 XMM-Newton fluxes are available (at similar or- 
bital phases), point to a quite constant X-ray flux versus signifi- 
cant Ha EW variations. If confirmed, this would disfavor a wind 
accretion scenario to explain the X-ray emission in the source. 
The constancy of the Nh points either to an emitter located at 
^ 10 12 cm from the compact object or to a stellar wind densities 
between 3 to 27 times lower than those expected from a spher- 
ically symmetric stellar wind model, since ionization does not 
seem to affect significantly the wind absorption. An upper-limit 
on the intrinsic A^h is also derived, being between 3 to 12 times 
smaller the one predicted by the mentioned stellar wind model at 
periastron passage. We propose that the multiwavelength emitter 
is not located inside the binary system, probably being the jet- 
like structure observed at radio frequencies. In any case, simulta- 
neous X-ray/optical observations are still required to clarify the 
stellar wind/X-ray link, and further XMM-Newton or Chandra 
X-ray observations are also necessary to further investigate the 
flux, photon index and intrinsic Nh evolution along the orbit. 
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